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We derive and evaluate expressions for the frictional Coulomb drag between dis- 
ordered two-dimensional electron gas layers. Our derivation is based on the memory- 
function formalism and the expression for the drag reduces to previously known 
results in the ballistic limit. We find that Coulomb drag is appreciably enhanced by 
disorder at low temperatures when the mean-free-path within a layer is comparable 
to or shorter than the layer separation. In high mobility two-dimensional electron gas 
systems, where the drag has been studied experimentally, the effect of disorder on 
the drag is negligible at attainable temperatures. We predict that an enhancement 
due to disorder and a crossover in the temperature-dependence and layer-separation 
dependence will be observable at low temperatures in moderate and low mobility 
samples. 
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I. INTRODUCTION 



Double-layer two dimensional electron gas (2DEG) systems, where electrons are confined 
to nearby parallel planes, are expected to exhibit many novel phenomena due to interlayer 
electron-electron interaction. For example, in strong magnetic fields Coulomb effects are 
expected to produce new incompressible ground states that exhibit the fractional quantum 
Hall effect |]T[] and to cause the collapse of certain integer quantum Hall effect gaps [Q. 
In zero magnetic field, it has been suggested that Wigner crystallization in double layer 
systems is favored by interlayer Coulomb interactions |^ and that in the case of electron- 
hole systems, excitonic superfiuity could result 0] from interlayer interactions. In recent 
experiments by Gramila et al. on electron-electron double layer systems [^,^ and in similar 
experiment by Sivan et al. on electron-hole systems [0 , the strength of interlayer interactions 
was studied directly by measuring the frictional drag of one two-dimensional electron gas 
layer on another. In these experiments a current fiowing in one layer tends to induce a 
current in nearby layers. If no current is allowed to fiow in the nearby layer an electric field 
develops whose infiuence cancels the frictional force between the layers. The transresistance, 
defined as the ratio of the induced voltage in the second layer to the applied current in the 
first layer, directly measures the rate at which momentum is transferred from the current 
carrying 2DEG to its neighbor. 

Drag between spatially separated electron systems due to Coulomb interactions be- 
tween carriers was first considered by Pogrebinskii and Price [^j. The experiments of 
Gramila et al. ||^,|^ and of Sivan et al. have stimulated recent theoretical attention, es- 
pecially to the case of drag between two-dimensional electron layers at low temperatures 
|^J^|l^|Tl||l^|l^JT^ . The interlayer Coulomb drag is caused by fiuctuations in the density 
of electrons in each layer since 2D layers with charge uniformly distributed will not exert 
any frictional forces upon each other. In this paper we examine for the first time the pos- 
sibility of enhanced frictional drag between disordered layers due to the diffusive nature 



of long- wavelength long-time electron density fluctuations. Disorder is known |15,Tq| to 
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enhance interaction effects and to lead to violations of Fermi liquid theory for individual 
two-dimensional electron gas layers. In the diffusive regime, where wavelengths are longer 
than the mean free path of the electrons and times are longer than the electron scatter- 



ing time, the electron density- density response function possesses a diffusion pole in 
momentum-frequency space. In perturbation theory the diffusion pole arises from dressing 
the electron-electron vertex [|l^] with corrections arising from impurity potential scattering. 



At shorter distances or shorter times the disorder vertex corrections are not important. We 
find that disorder enhances the interlayer drag at low temperatures, changing the tempera- 
ture dependence of the drag from to — InT. In very high-mobility samples (on which 
existing experiments have been performed) or in samples with small layer separations the 
corrections due to disorder become important only at extremely low temperatures. For 
samples with lower mobility or more widely separated layers the influence on the interlayer 
scattering rate from disorder scattering should be easily measurable. 

Most [^|TO],|rT| previous work on interlayer friction has been based on Boltzmann trans- 



port theory which cannot capture disorder enhanced interaction effects. (An exception is the 
work of Vasilopoulos and co-workers [0].) In Section II we present a derivation of the expres- 
sion for the frictional transresistance based on the memory function formalism which is 



sufficiently general to treat the case of disordered 2DEG layers. Although phonon-mediated 
P JI^JTOd interlayer interactions are not discussed explicitly in this paper the expression we 
derive in Section II are sufficiently general that other coupling mechanisms can be incorpo- 
rated as an effective-interlayer interaction potential. Similarly particle-particle interaction 
vertex corrections, which are probably quantitatively important especially for electron-hole 
double layer systems 0, can also be incorporated into the results we derive in Section II 
as a contribution to the effective electron-electron interaction |jl8[. Readers interested only 
in the application to disordered double-layer systems should proceed to Section III where 
we discuss how disorder within the layers influences the interlayer friction. The interlayer 
scattering rate of typical electron-electron double layer samples is evaluated numerically in 
Section IV. We show that the temperatures below which disorder becomes important de- 



creases very rapidly with increasing mobility. A brief summary of our findings concludes the 
paper in Section V. 



II. MEMORY FUNCTION FORMALISM DERIVATION 

The memory function formalism provide a very convenient method for deriving a flexible 
expression for the transresistance of double layer 2DEG systems. The Kubo current-current 
correlation function formula for the conductivity is first converted into a force-force correla- 



tion function expression for the resistivity by making use of Mori's projection operator. 
This force-force correlation function is then evaluated at lowest order in the screened inter- 
layer interaction to obtain an approximate expression for the transresistance. The advantage 
of using the force-force correlation function rather than the current-current correlation func- 



tion is that it yields a reasonable approximation even when evaluated at lowest order pl| , |22 
At this level the results are physically equivalent to momentum-balance |JI^ approximations 
or to relaxation-time approximations in a Boltzmann-transport approach The deriva- 
tion for the situation of present interest is sketched in the following paragraphs. 

For notational simplicity we restrict ourself to the case of zero magnetic field so that the 
current is in the same direction as the applied electric field and the conductance, therefore, 
forms a 2x2 matrix with respect to the layer indices. (The final expression for the transre- 



sistance is equally valid p4[ in the presence of a magnetic field.) We need to consider only 
the long wavelength limit of the conductance. To use the memory function formalism it is 
convenient to write the Kubo formula in the form 

a,,{u) = ^j^ dte"^\Ut)Jj) (1) 

where (3 = l/ksT, v is the cross section area of the 2DEG layers, and the indices % and j 
are layer labels. OijiuS) gives the current density induced in layer % due to an electric field 
in layer j. J is the zero wavevector Fourier component of the current density operator. The 
inner product appearing in Eq. is defined by 



4 



CAB{t) = {A{t),B) 



p-^ f'^ d\{A\t),B{ih\)) (2) 







In Eq. (H) the angle brackets denote thermal averages. The following relationship can be 
used to change Eq. (Q) into the more familiar form of the Kubo formula 



mCAB{t) = ^{[A{t),B]). (3) 

The projection operator method is now used to obtain an expression for the matrix 
inverse in layer indices of Cj-jy We define a superoperator V which 'projects' an operator 
O onto the current density, and its complement Q, by 

p6=$:M^ (4) 

k {Jk, Jk) 

vd = d-Qd (5) 

It is useful to define a matrix Xij'- 

X., = ^C^J.J,(0) = (4i,) = ^5.,. (6) 
Here is the areal density of 2D electrons in the z-th layer. Following the usual development 



of the memory function formalism |13 , we obtain an equation for the matrix inverse of the 
Fourier transform of Cj^j.{t): 

[CJ^J^{z)]-' = ^X-\-^zl + M{z)] (7) 

where 



A dot over an operators denotes its time derivative at t = and the Liouville superoperator 
is defined by CO = [H,0], where H is the Hamiltonian. In obtaining Eq. (|^), the time 
reversal invariance condition VJ = has been applied. Combining Eq. (|I|), Eq. (|^ and 
Eq. (H) we obtain an expression for the resistivity matrix: 



■'■^ V Jo 



V Jo 

Pij relates the electric field in layer i to the current density in layer j. With the relation Jj 
—e/mFi, we obtain the force-force correlation function expression for the trans-resistance 



poo 

/ dte''\FL,e-'Q^'Fn). (10) 
Jo 



PLR — r, 

nL'TiRe^i' Jo 

It is easy to demonstrate that plr is identically zero in the absence of interlayer coupling, 
since the forces in left and right layers are uncorrelated. The leading contribution to the force 
operator from interlayer interactions can be expressed in terms of the interlayer interaction 
potential Ue{q) and the electron density g{q) 



^RiL) = ±-Y.^QL{q)QR{-q)Ue{q). (11) 



To leading order in interlayer interactions e~*^^*Fi in Eq. (|10|) can be replaced by e^*^*Fj. 
This replacement leads to the desired force-force correlation function expression for the 
trans- resistance 

Plr{z) = ^ r dte'^\FL{t), Fn)o. (12) 

nLnRe'^u Jo 

The subscript on the inner product in Eq. ([T2|) indicates that it should be evaluated in 
the absence of interlayer interactions. Substituting the explicit expression for the interlayer 
forces gives 



with 

Mt) = QLi-q,t)gRiq,t). (14) 



The correlation function for decoupled layers appearing in Eq. (|T3]) is related to the 
isolated layer density fluctuations. Using a representation of exact eigenstates that in the 
z = limit it is easy to show that 
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" 2n^^Je2 / dujSL{q,uj)SR{-q,-uj) (15) 

where Si{q,uj) is the dynamic structure factor for layer i: 

S,{q,u;) = -Y.exp{-(3En)\{n\p,{q)\m)\^6{u;-{E^-En)/n) (16) 

^ n,m 

It is usually more convenient to express the resistance in terms of individual layer response 
functions rather than the dynamic structure factor. We relate Si{q,u!) to the retarded 
density-density response function for layer i, XiiQ, ^) by applying the fluctuation-dissipation 
theorem [|l^ 



This gives us the final form of our expression for the drag resistivity, which is summarized 
diagrammatically in Fig. (^. 



Plr 



irriLnRe^ ^ Jo ^^hu; ^ ^-pn^ _ 2 ^ ^ 



III. DISORDER AND SCREENING 

The electron density-density response function xiQ,^) in Eq. (|18|) can be obtained by 
applying many-body perturbation theory methods to a 2DEG whose Hamiltonian contains 
disorder and (or) interaction terms. For a non-interacting disorder-free 2DEG the response 
function can be evaluated analytically p5[. Disorder leads to an enhancement in lmx{q,uj) 



at low frequencies and long wavelengths and, as we discuss below, can enhance the interlayer 
friction. The enhancement reflects the increased spatial correlation of states with nearby en- 
ergies in disordered systems. The effect of disorder on the friction can be described without 
making a specific model of disorder by invoking the Einstein relation between the conductiv- 
ity and the diffusive density-density response at long wavelengths and low frequencies ^6 



We introduce a phenomenological intralayer electron (transport) scattering time r. r is re- 
lated to the mobility by /i = er/m and at low temperatures is related to the mean-free-path 



by / = rhkf/m. For gZ > 1 or ur > 1 we assume that disorder is unimportant and approx- 
imate the density-density response function by the non-interacting electron resuh [^]. At 
zero temperature and zero disorder 

= ^^{- - C^mq/mf -{C0 + e,f/nr' - C-\{kfq/mf - {u - e,f 1%^'} (19) 

where Eg = hq'^/2m, /i is the chemical potential and the C± are: 

C± = sign(£q ± cj) if {kfq/mY - {uj±ey/h'^ <0, 
C± = ±i if {kfq/mf - (cj ± Eqf/h'^. > 

However for ql < 1 and ujt < 1 disorder becomes important. In this diffusive regime 
the electron density- density response function is completely characterized by the diffusion 
constant [1^: 



x^iq,^) = xFiq,^) 

dn Dq^ 



(20) 



(i/x Dq^ — iuj 

q <l/l uj <1/t 

where D = 1^/2t is the diffusion constant. 

The derivation of the expression for the transresistivity in the previous section is valid up 
to second order in the interlayer interaction. For the system of physical interest the interlayer 
interaction is Coulombic and it is essential to include screening in order to get qualitatively 
correct results. In this paper we adopt the usual expediency of employing the second- 
order expression with the interlayer interaction replaced by a screened interlayer interaction 
and argue that this includes the most important higher-order effects. In the disorder free 
limit our expression for the transresistivity then becomes identical to those derived using 
other approaches in earlier work [|5|, p!0| , pJ] , p^ , p^ . The random-phase-approximation (RPA) 
screened interlayer interaction is 
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^ [l + K(g)xL(g,c.)][l + K(g)x/?(9,^)]-v;'(g)xL(g,u;)x«(g,^) ^ ^ 

where the bare intra- and inter-layer electron-electron interaction potentials are Va{q) = 
27re^/g and Ve{q) = Va{q)e~'^'^ where d is the separation between the layers. In the above 
expression either the ballistic or the diffusive form for Xl{r) should be used as appropriate. 
Note that the interlayer interaction is cut-off by the factor e~'^'^ for q > 1/d. The layer 
separation dependence of the friction in both diffusive and ballistic limits results from this 
cutoff. Physically the cutoff reflects the fact that charge fluctuations in one layer with a 
wavelength shorter than the layer separation get averaged out when viewed from the other 
layer. 

The expression for the trans-resistance of Eq. (|T^) can be split into contributions from 
the ballistic and diffusive regimes. With p^^ = nue^TLji/m, we have 

TrI = + r^' (22) 



where 



2TT'^mnL Jo Jo e'^^'^ + e"^^"^ - 2 



and 



r.^ - Mi t' Mu /^^^i^^ (24) 



with 



2'K'^mnLJo Jo ^(Shu ^ ^-fsn^ _ 2 



Tj^^ is the result for a disorder-free 2DEGs and r^^ is the correction due to the the enhanced 
fluctuations at long wavelengths and low frequencies in disordered systems. In the next 
section we discus the evaluation of these expressions. 

IV. NUMERICAL RESULTS AND DISCUSSION 

The dependence of the interlayer scattering rate on temperature and on layer separation 
depends on whether the interlayer scattering is dominated by r^^ or r^^^. In Fig. (§) and 
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Fig. we show numerical results for tb and functions of temperature calculated 

for two different values of layer separation for a high mobility two-dimensional electron gas 
sample. The data in these figures are obtained from numerical evaluation of Eq. (^) and 
Eq. (^) with the input parameters taken from the experiment of ref. 0. In Fig. (^) we see 
that r^^ ~ T'^/d^ at low temperatures, as pointed out in Ref. ||^. (The d~'^ dependence can 
be recognized by noticing that the scattering rate decreases by a factor of approximately 
four when the layer separation increases by a factor of v^-) From Fig. @ one can see that 
falls off more slowly with both temperature and layer separation as T ^ 0. The inset 
to Fig. (^) establishes that for T <^ Tr, ~ — T^lnT/rf^. The dimensionless temperature 
scale for Fig. (^) is the Fermi temperature (Tp = Ep/ks) which is about 60K for this 
sample, while the dimensionless temperature scale for Fig. (|]) is the disorder temperature 
{Tt- = h/ksr) which is ~ 56mK for this sample. Comparing Fig. (0) and Fig. (^) we see 
that for this high-mobility sample the disorder correction is smaller than one part in 10^ at 
temperatures above ~ lOmK. 

The origin of the temperature and layer dependence seen in Fig. and Fig. can be 
understood by looking at the limit of large layer separations where the interlayer scattering 
rates can be evaluated analytically. The evaluation of t^^ at large layer separations and 
low temperatures for the disorder-free limit, where the relation tb oc T'^/d^ holds, has 
been carried out previously by several authors P, p!0| , pT| . (The results quoted in Ref. ( 



0) and Ref. ( [|1^]) are in error by a factor of two.) We rederive those results here to 
allow a comparison with the disordered case. For T <^ Tp and d ^ kp^ only the low 
frequency and long wavelength limit of Imx contributes importantly to Eq. ([T8|). {kp is 
the Fermi wavelength.) From Eq. (|T9|) it follows that in this limit liexiq,^^) — dn/dfi, and 
Imx(g, cij) = dn/dfi{2huj/Ep){kp/q). For T <^ Tp, we can ignore the contribution of Imx 
to screening the interlayer interaction and it follows from Eq. (pT]) that we can replace the 
interlayer interaction by 

^'^^^^ kj,psmh.qd ^ 
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where kTF = 27re'^ dn/dfi is the single layer Thomas- Fermi screening wavevector. Note 
that for q ^ the effective screening wavevector is 2k'^pd , which is proportional to the 
layer separation. (For GaAs kxF = 0.2nm~^ independent of electron density.) With these 
approximations the integral over frequency and wavevector are known and we obtain 

-7rC(3)(A;BT)2 



(27) 



^ lQhef{kTFdY{kFdy' 
In this result two powers of d~^ may be associated with the enhanced screening of the 
interlayer interaction at large separations and two powers of d~^ with the combination of 
phase space considerations which cause the integrand to to vary as for small q. For layer 
separations larger than the mean free path the electron response is diffusive over the entire 
range of wavevectors contributing importantly to Eq. (|l^). For low temperatures it follows 
from Eq. (pO|) that we may use Rex(q', tu) = dn/dfi and Imx(g, cj) = dn/ dii[ujDq^ / {oj'^ + 
[Dq^Y)]. Again we may ignore the contribution to screening from Imx so that the relevant 
limit of the screened interlayer interaction is unchanged. At small tu, Imx oc compared 
to the q~^ of the ballistic case. The integrand of the wavevector integral thus goes as q^^ at 

1 /2 

small wavevector. This logarithmically divergent wavevector integral is cutoff at q ^ u/D ' . 

The remaining frequency integral is elementary and we obtain 

, ^ -n{kBT)MT/Tr) 
^ UhefiqTFdyikFiy ^ ^ 

The change in the layer separation from ci"^ in the ballistic case to d~^ in the diffusive case 

can be traced directly to the change in the wavevector dependence of Imx from q~^ to q~^. 

For d <^ / Eq. (|28|) correctly gives the contribution to the drag from q <^ l^^ and Eq. (pT]) 

gives the contribution from l^^ < q < d^^. Because of the different temperature dependence 

it is still true that the contribution from the diffusive regime will dominate at sufficiently low 

temperatures. Comparing Eq. ( pH]) and Eq. ( pTj) we can estimate the crossover temperature: 



T,~T.exp[-3(//rf)V4C(3)] (29) 

For high-mobility samples the diffusive enhancement of the drag will be observable only at 
extremely low temperatures. For example, in GaAs is about 0.2K and / ~ lOyU for a 
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sample with a mobility of ~ lO^cm^ / sV and a typical density. For a layer separation of 
~ 500^4° this implies that ~ 10^^^^ K. For the samples in the experiment of ref. 0, 
which are of extremely high mobility and small layer separations, the value of is 56mK 
and the d/l is about 10~^. r^^ is smaller than r^^ by a factor of 10^'"^ at T ~ T^. The 
correction term, r^^, will be difficult to observe for accessible temperatures. In samples with 
lower mobility and/or thicker barriers between the layers, r^^ and t^^ have comparable 
amplitudes for T ^ Tr and the contribution from the diffusive regime dominates at lower 
temperatures. In Fig. we plot the relative contribution to the drag from the diffusive 
regime vs. mobility for a layer separation of 50nm and n = 1.5 x lO^^cm^^. These results 
show that the effect of disorder will become easily observable at typical low temperatures 
for samples with mobilities below ~ lO^cm^V^^s"^. 

It is possible to fabricate double-layer systems in which one layer is much more disordered 
than the other. In particular one may have / ^ d for one layer and / < d for the other layer. 



Following the same steps leading to Eq. (|28|), it is possible to derive an expression for the 
low-temperature transresistance in such a system by using one diffusive response function 
and one free electron response function: 



Tt 



(30) 



12hef{kTFdY{kFd){kFiy 
The dependences on temperature and layer separation, r^"^ ~ T"^ /d^, are easily understood 



by comparing it to Eq. p7| ) and Eq. (^) and noticing that the integrand in this case 
approaches a constant at small wavevector transfers. In Fig. (^, results are shown for the 
relative correction due to disorder enhancement for the case where one layer consists of free 
electrons while the other layer has a finite mobility. The relative correction is essentially 
independent of temperatures at low temperatures and it is weaker than in the case where 
both layers are disordered. 
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V. SUMMARY 



Using the memory-function method, we have derived an expression for the transresis- 
tance of double layer systems which is sufficiently general to treat the case of disordered 
layers. The expression has been evaluated as a function of temperature, layer separation (d) 
and in-plane mobility /i. Both the case where only one layer is disordered and the case where 
both layers are disordered have been considered. We find that the drag varies as c^"^, 
and at low temperatures for clean, single-layer disorder, and double-layer disorder cases 
respectively. In the case of double-layer disorder the transresistance varies as — In(T') at 
low temperatures, otherwise the transresistance is proportional to T^. The low-temperature 
drag is proportional to and for double-layer disorder and single-layer disorder re- 
spectively. Except for the case of extremely low temperatures the crossover from clean to 
disordered regimes occurs when the mean-free-path within a layer becomes smaller than the 
layer separation. In very high mobility two-dimensional electron gas systems, where the 
transresistance has been studied experimentally up to the present, the effect of disorder on 
the drag is negligible at available temperatures. We predict that the enhancement due to 
disorder and the associated crossovers in the temperature-dependence and layer-separation 
dependence will be observable at low temperatures in moderate and low mobility samples. 
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FIGURES 

FIG. 1. Schematic diagrammatic rendering of the memory function expression for the transre- 
sistivity to lowest order in inter layer interactions. (Eq. (|l^)) 

FIG. 2. The interlayer scattering rate for a disorder-free 2DEG samples as a function of temper- 
ature, t = T/Tp where Tp is the Fermi temperature. The solid curve was obtained for d = 315^° 
and the dashed curve was obtained for d = \/2 x 315^°. These curves were calculated for typical 
and equal densities in the two layers: ul = ur = 1.5 x IQ^'^cm^^V so that Tp ~ QQK. 

FIG. 3. The correction to the interlayer scattering rate due to disorder-enhanced interactions 
at long wavelengths and low temperatures. Here t = T jT^-. The insert shows that oc T^lnT/d^ 
at low temperatures. The solid curve is for d = 315^1° and the dashed curve is for d = \pl x 315^°. 
Other parameters used are taken from the experiments of Ref.( [^): ni^ = n/j = 1.5 x lO^^cm"'^, 
and /X = 3.4 x \^^cm} j sV . The corresponding Tj- = 56mK° so that disorder corrections would not 
be observable at available temperatures. 

FIG. 4. Relative correction to the interlayer scattering rate due to disorder-enhanced inter- 
actions for several temperatures as a function of sample mobility. These results were obtained 
for a typical layer separation: d = 500^°. The dashed line is for T = l.OK, the dotted line for 
T = O.IK and the solid line for T = O.OIK. The density was taken from the experiments of 
Ref.( [|[). For these temperatures the diffusive contribution become important for mobilities below 
about 10^ cm^s ^. The marked lines show the result for the case where one layer has a finite 
mobility while the other layer has an infinitively large mobility. 
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